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ABSTRACT 

We  present  a  ful  1-bands tructurc  calculation  of  temperature-  and  wavelength-dependent 
two-photon  absorption  (TP A)  coefficients  and  free-carrier  absorption  (FCA)  cross 
sections  in  InAs,  InSb,  Hgu-Cd.Te  alloys.  Although  the  wavelength  dependence  of  the 
TP  A  coefficients  agrees  well  with  widely  used  analytical  expression,  our  calculated 
values  are  smaller  by  a  factor  of  1.2  to  2.5.  In  addition,  the  TP  A  coefficient  depends 
sensitively  on  the  photoexcitcd  carrier  density  in  small  gap  materials.  The  FCA  is  found 
to  be  due  predominantly  to  holes.  The  FCA  cross  section  is  independent  of  the  carrier 
density,  but  strongly  dependent  on  the  temperature.  The  calculated  coefficients  and 
lifetimes  are  fitted  to  closed-form  expressions  and  used  in  solving  the  rate  equations  to 
obtain  the  transmitted  pump  and  probe  intensities  as  functions  of  incident  intensity  and 
sample  thickness.  The  calculated  pump  transmission  and  time-dependent  probe 
transmission  in  InAs  agree  very  well  with  the  measured  values. 
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INTRODUCTION 


The  nonlinear  optical  absorption  and  refraction  properties  affect  the  propagation 
of  high  intensity  light  through  transparent  materials.1'2  Several  physical  mechanisms 
including  two-photon  absorption,  free-  carrier  absorption,  thermal  and  free  carrier 
refraction  are  responsible  for  nonlinear  properties.3'6  Both  theory  and  experiments  are 
constantly  being  improved  to  accurately  evaluate  nonlinear  absorption  and  refraction 
coefficients.7"22  When  the  energy  of  high-intensity  incident  photons  is  more  than  a  half  of 
the  band  gap  of  the  semiconductor,  the  electrons  in  the  valence  band  (VB)  absorb  two 
photons  to  reach  unoccupied  states  in  the  conduction  band  (CB).  Then  the  photoexcited 
electrons  and  holes  further  absorb  a  photon  with  or  without  assistance  of  a  phonon, 
causing  further  reduction  in  intensity  as  the  light  travels  through  the  material.  Since  the 
two-  photon  absorption  (TP A)  plays  a  vital  role  to  initiate  on  the  nonlinearity  in  the 
material,  a  number  of  studies  have  been  undertaken  to  accurately  evaluate  the  TPA 
coefficient,  However,  the  extraction  of  TPA  values  from  experiments  often  depended  on 
the  ability  to  control  other  factors  that  simultaneously  affect  the  nonlinear  absorption. 
Consequently,  the  extracted  values  varied  substantially.7" 10,1 5,1 8-21  On  the  theoretical  side, 
the  Wherrett’s  expression12  with  nonparabolic  correction11  is  used  widely.  Later  it  was 
modified  and  the  interaction  matrix  element  was  parameterized  based  on  the  k.p  models 
to  successfully  predict  the  trends  in  various  semiconductors.14"1  This  expression  has 
served  an  excellent  role  in  providing  qualitative  trends  and  giving  approximate  values  for 
the  TPA  coefficients  in  understanding  light  propagation.  However,  an  accurate 
evaluation,  with  full  band  structures  and  matrix  elements  included,  of  TPA  coefficient 
made  recently22,23  has  revealed  that  the  calculated  values  differ  substantially  from  those 
given  by  the  simple  analytical  formula. 

Light  propagation  is  further  affected  by  free  carrier  absorption  (FCA),  in  which 
the  TPA  generated  photoexcited  carriers  absorb  photons.  The  cross  section  of  the  free 
carrier  absorption,  defined  as  the  ratio  of  free  carrier  absorption  coefficient  to  the  free 
carrier  density,  is  in  general  assumed  to  be  a  constant  in  light  propagation  studies/"6 
However,  in  our  detailed  calculations  we  find  that  the  cross  section  is  strongly  dependent 
on  wavelength  and  temperature.  We  have  carried  out  these  calculations  with  full  band 
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structures  for  InSb,  InAs,  and  Hgi,xCdxTe  alloys.  The  room -temperature  band  gap  of 
InAs  and  InSb  are  0.35  eV  and  0.175  respectively. 

In  this  paper,  we  first  provide  a  brief  description  of  the  bandstructures,  the  quasi 
Fermi  energies  and  nonlinear  absorption  coefficients  (TPA  and  FCA).  The  wavelength 
and  carrier  density  dependence  of  the  TPA  and  FCA  coefficients  in  InAs,  InSb,  and 
HgCdTe  alloys  are  reported.  The  calculated  values  of  TPA  and  FCA  in  InSb,  InAs,  and 
HgCdTe  alloys  (with  x=0.238  and  x=0.344)  are  also  fitted  to  simple  analytical  form  so  it 
can  be  used  easily  in  light  propagation  modeling.  3‘6’24  These  values  are  used  in  obtaining 
a  solution  to  rate  equation  that  describes  the  light  propagation  through  InAs.  The  results 
are  compared  with  experiments. 

THEORY 

Second-order  perturbation  theory  is  used  to  calculate  the  TPA  and  FCA 
coefficients22  and  Auger  and  radiative  recombination  lifetimes.2'’  The  dipole  matrix 
elements,  required  to  study  light  interaction  with  matter,  and  the  scattering  matrix 
elements,  required  in  lifetime  calculations,  are  calculated  from  underlying  bandstructures. 
Care  is  taken  to  describe  the  occupancy  of  initial  and  final  states  with  appropriate  Fermi- 
Dirac  distributions  and  to  include  accurate  sum  over  crystal  wave  vectors  spanning  the 
entire  Brillouin  zone  (BZ). 

Band  Structures 

The  calculations  reported  here  are  based  on  a  hybrid  pseudopotential  tight-binding 

26 

Hamiltonian.  In  this  formalism,  the  tight-binding  Hamiltonian  contains  two  parts.  The 
dominant  part  is  derived  from  the  empirical  pseudo  potential  and  includes  all  long-range 
interactions.  A  short-range  empirical  Hamiltonian  is  added  to  fine-tune  the  calculated 
band  structure  to  agree  with  the  measured  energy  gaps  and  effective  masses  at  various 
symmetry  points  in  the  BZ,  In  the  spJ  basis,  with  the  spin-orbit  coupling  included,  the 
Hamiltonian  takes  16x16  matrix  form.  In  the  case  of  alloys,  we  find  that  the  correction  to 
bands  from  the  alloy  disorder  is  small  near  the  band  edge,  and  hence  the  virtual  crystal 
approximation  can  be  used  to  obtain  the  alloy  band  structures.  The  resulting  Hamiltonian 
matrix  is  diagonalized  to  yield  the  eigenvalues  (energies)  and  eigenvectors  (wave 
functions). 
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Since  a  high  intensity  light  absorbed  in  a  material  can  raise  the  material  temperature 
substantially  above  room  temperature,  we  are  interested  in  using  bandstructures  over  a 
wide  range  of  temperatures  from  77  K  up  to  -500K..  Although  the  temperature  (T) 
dependence  of  the  matrix  elements  can  be  included^'  for  the  calculation  of  the  band 
structure,  those  calculations  are  extremely  time-consuming,  particularly  when  the 
integration  over  the  full  BZ  is  required  in  lifetime  calculations.  Here,  as  an  alternative, 
the  effect  of  T  on  the  Hgi^Cd/Te  alloy  band  structure  is  included  indirectly  as  follows. 
First,  for  a  given  T  and  Cd  concentration  x ,  the  band  gap  Eg  in  eV  is  determined  from  the 
widely  used  relationship,28  -0.302  +  1.93jc  +5.35  10'4  (l-2x)*T-0.81x2+0,832x3.  We  then 
adjust  the  value  of  x  in  our  zero-temperature  Hamiltonian  to  x'  to  get  the  band  gap  of  Eg 
The  zero-temperature  Hamiltonian  with  the  value  of  x'  is  then  used  in  further 
calculations.  In  the  case  of  InAs29  and  InSb30  compounds,  the  parameters  of  the  short- 
range  empirical  Hamiltonian  are  adjusted  at  each  temperature  to  yield  the  band  gap  in  eV 
given  by  0.415  -  2.76x1  O'4  T2  /(T+83)  and  0.24  -  6x1  O'4  T2/(T+500)  respectively.  The 
calculated  band  gap  energy  as  a  function  of  T  is  shown  by  dashed  line  in  Figure  1  (a) 
InAs,  (b)  Hgo.66Cdo.34Te,  (c)  InSb,  and  (d)  Hgo.76Cdo.24Te.  We  see  that  the  band  gap 
decreases  with  T  in  InAs  and  InSb,  but  increases  in  HgCdTe  alloys. 

Fermi  Level 

The  study  of  electron-photon  interaction  requires  the  knowledge  of  the  electron 
density  and  its  distribution.  The  Boltzmann  distribution  is  often  used.  However,  at  high 
temperatures  and  for  small  gap  materials,  the  bands  are  nearly  or  fully  degenerate  and  the 
Fermi-Dirac  (FD)  distribution  is  required  in  the  functional  form, 

f(E,Ef)  =  [l  +  e(E"Ef)/kBT]"' 

where  kB  is  the  Boltzmann  constant.  The  equilibrium  Fermi  level  (FL)  or  Fermi  energy  Ef 
is  obtained  from  the  charge  neutrality  condition  that  the  total  positive  charges  (the  sum  of 
hole  density  p0  and  ionized  donors  Np )  equals  the  total  negative  charges  ( the  sum  of 
electron  density  no  and  ionized  acceptors  N“  ).31  The  quantities  no,  po,  N^and  NA  in 
turn  depend  on  the  Fermi  energy: 
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n0=2  Jpc(E-Ec)f(E,Ef)dE;  p0  =  2  Jpv(E-Ev)f(E,Ef)dE; 


Under  high-intensity  of  light  illumination,  very  large  photoexcited  electron  density 
(An)  and  hole  density  (Ap)  are  created.  Although  An=  Ap,  they  can  be  far  above  their 
intrinsic  values  of  no  and  po-  Consequently,  the  TPA,  FCA  coefficients  and  lifetime 
calculations  have  to  be  carried  out  not  only  under  degenerate  conditions,  but  also  under 
nonequilibrium  condition  in  which  the  electrons  in  the  CB  and  the  holes  in  the  VB  have 
different  Fermi  energies.  These  quasi  Fermi  energy  for  electrons  (Ef )  and  holes  (Ef )  are 
determined  from  the  following  conditions: 


n0  +  An  -  2  Jpc(E-  Ec)f(E,Ef  )dE 


Ev 


p0+Ap  =  2  Jpc(E-Ev)f(E,Efv)dE 


where  Ee  and  Eu  are  the  band  energies  at  CB  minimum  and  and  VB  maximum 
respectively.. 

In  Figure  1,  we  have  shown  the  calculated  electron  (thin  solid  line)  and  hole  (thick 
solid  line)  Fermi  energies  as  functions  of  the  lattice  temperature  and  excited  carrier 
density  in  (a)  InAs,  (b)  Hgo.66Cdo.34Te,  (c)  InSb,  and  (d)  Hgo.76Cdo.24Te.  When  the 
excited  electron  (An)  and  hole  (Ap)  densities  are  zero,  we  obtain  the  intrinsic  Fermi 
energy,  shown  by  dotted  line.  All  energies  are  measured  with  respect  to  the  top  of  the  VB 
at  that  temperature.  At  low  temperature,  the  intrinsic  FL  is  near  mid-gap  and  moves 
closer  to  the  CB  edge  (shown  by  dashed  line)  with  an  increase  in  T.  In  the  case  of  InSb, 
we  see  that  the  FL  lies  above  the  bottom  of  the  CB  at  high  temperature,  making  the 
material  fully  degenerate.  In  HgCdTe  alloys,  the  band  gap  increases  with  T  and  hence  the 
intrinsic  FL  is  lies  below  the  CB  for  all  the  temperatures  considered  here.  When  the 
photoexcited  carrier  density  (An=Ap)  is  finite  and  larger  than  the  intrinsic  carrier  density 
at  that  T,  the  quasi  FL  for  electron  and  holes  are  closer  to  the  CB  and  the  VB 
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respectively.  The  quasi  FLs  merge  with  the  intrinsic  FL  when  the  intrinsic  carrier  density 
is  larger  than  the  photoexcited  carrier  density.  The  intrinsic  carrier  density  increases  with 
T  and  decreases  with  Eg,  as  seen  in  Figure  2.  Consequently,  quasi  FLs  play  major  role 
even  in  systems  with  moderate  excited  densities  at  low  T  and  in  systems  with  large 
densities  at  high  T  in  general.  This  trend  is  clearly  seen  in  Figure  1 .  Recently,  the  variable 
Hall  studies,  in  combination  with  multiple  carrier  analysis,  have  demonstrated  to  be  able 
to  a  yield  accurate  measurement  of  the  carrier  density  in  doped  samples.'*2'34  The  carrier 
densities  calculated  using  our  band  structures  were  found  to  be  in  excellent  agreement 
with  the  measurements.17 

Two-photon  absorption 

Electrons  in  the  heavy-hole  (HH)  and  light- hole  (LH)  valence  bands  absorb  two 
photons  to  reach  the  CB  states.  In  the  second  order  perturbation  theory,  the  absorption 
rate  is  proportional  to  the  occupation  of  the  initial  state  in  the  VB,  availability  of  the  final 
state  in  the  CB,  and  the  dipole  matrix  element  between  the  VB  and  CB  states  at  a  given 
crystal  momentum.  Explicitly,  the  expression  for  TPA  coefficient  p  is: 


P  = 


1  ( 27ieY 

TIC  y  m  J 


1 

(®0n)3 


£(i-f(E;,E?))|i<f 


f(E‘,E?)s(E‘-Ej +:»<»,,) 


k  _  y  y  <4>p  I  VP  !<!>?■  )(<(>?■  Iej  -pM>a) 

“P  ~Ts  e£-e;.+jk 

where  Ejj  (<jip )  is  the  electron  energy  (wave  function)  in  the  band  p  at  the  wave  vector  k 

in  the  BZ,  e  and  m  are  free  electron  change  and  mass,  n  is  the  real  part  of  the  refractive 
index,  p  is  electric  dipole  vector,  ej  is  the  polarization  of  light,  and  fte>0  is  photon  energy. 

Note  that  the  FD  distribution  for  the  electrons,  f(Ep )  uses  appropriate  quasi  Fermi 


energy  Ep  and  E)  depending  on  whether  p  denotes  a  valence  or  a  CB.  In  the  above 
expression,  a  and  p  denote  respectively  initial  and  final  state  in  any  occupied  VB  and  any 
unoccupied  CB  that  satisfy  the  energy  conservation  condition  specified  by  the  8  function. 
A  careful  analysis  shows  that  the  intermediate  state  p  can  be  any  of  the  unoccupied  CB 
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state  and  one  VB  state  from  which  the  electron  excitation  begins.  The  sum  over  k  is 
considered  over  entire  the  BZ. 

We  used  full  bands  to  accurately  accurate  FD  distributions,  band  energy,  wave 
functions,  and  optical  matrix  elements,  to  obtain  p  as  a  function  of  the  photon  wavelength 
and  T.  The  calculated  p  at  T=300K  and  77K.  are  shown  in  Figure  3  for  (a)  In  As,  (b) 
Hgo.66Cdo.34Te,  (c)  InSb,  and  (d)  Hgo.76Cdo.24Te.  When  the  wavelength  is  longer  than 
twice  the  cutoff  wavelength,  electrons  in  the  VB  cannot  reach  the  CB  and  p  drops  to  zero. 
As  the  wavelength  is  decreased,  more  states  are  available  and  p  increases.  However,  with 
further  decrease  of  the  wavelength,  the  energy  difference  between  the  initial  and  Final 
states  increases  and  the  rate  calculated  by  the  second-order  perturbation  theory  decreases, 
resulting  in  a  decreased  p.  This  trend  is  seen  in  all  materials  considered  here.  As  the  one 
photon  absorption  will  start  to  dominate  at  wavelengths  shorter  than  the  cutoff 
wavelength  for  the  given  T,  the  values  of  p  calculated  in  that  spectral  region  is  not  useful. 
In  InAs  and  InSb,  the  band  gap  decreases  with  T,  and  consequently  p  increases.  For 
similar  reasons,  p  in  HgCdTe  decreases  with  T.  Although  the  band  gaps  of  InAs  and 
Hgo.66Cdo.34Tc  at  room  temperature  are  equal,  p  in  HgCdTe  alloy  is  considerably  higher 
than  in  InAs,  owing  to  the  difference  in  the  density  of  states. 

For  an  intrinsic  material  at  a  very  low  temperature,  the  initial  state  is  full  and  the 
final  state  is  empty.  However,  the  occupancy  or  availability  of  the  states  depends  on  T 
and  the  excess  carrier  density,  given  by  FD  distribution  as  discussed  in  the  previous 
section.  Consequently,  the  calculated  p  has  a  dependence  of  the  photoexcited  carrier 
density  as  shown  in  Figure  4  at  (a)  300IC  and  (b)  77K.  At  300K,  the  intrinsic  carrier 
density  in  the  large  band  gap  InAs  and  Hgo.66Cdo.34Te  is  small  and  the  calculated  p  shows 
only  slight  variation  with  photoexcited  carrier  density.  However,  in  the  small  gap  InSb 
and  Hgo.76Cdo.24Te,  the  intrinsic  carrier  density  is  large.  Additional  photoexcited  carriers 
deplete  (fill)  the  near  band  edge  VB  (CB)  states.  Hence,  the  occupancy  (availability)  of 
initial  (final)  state  is  different  from  one  (zero),  resulting  in  a  decreased  p.  We  see  that  the 
calculated  p  start  to  decrease  for  the  photoexcited  carrier  densities  in  excess  of  mid  1016 
cm"1.  This  effect,  commonly  known  as  Moss-Burstein  effect,  is  dominant  for  small  band 
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gap  materials  and  at  low  T.  The  |3  calculated  at  77K  shows  even  stronger  carrier  density 
dependence,  as  shown  in  Figure  4b. 

It  is  instructive  to  compare  the  calculated  p  to  that  obtained  by  most  commonly  used, 
but  approximate,  methods.  Quantitatively,  the  calculated  wavelength  dependence  agrees 
well  with  that  predicted  by  the  Wherrett’s  or  Van  S  try  land’s  expression.8,11  However,  our 
values  calculated,  at  4.8  pm  and  low  carrier  density,  are  smaller  than  those  obtained  with 
Van  Stryland’s  expression  by  a  factor  of  2.5  for  Hgo.66Cdo.34Te  and  1.5  for  InAs. 
Similarly,  the  values  calculated  at  9.6  pm  and  low  carrier  density  are  smaller  by  a  factor 
of  2.2  for  InSb  and  1.2  for  Hgo  74Cdo.2fiTe.  Although  the  values  calculated  by  either 
method  are  within  the  experimental  accuracy,  it  is  important  to  note  that  the  newly 
identified  N-  and  T-  dependence  will  be  important  in  light  propagation  studies. 

The  wavelength-dependent  calculations  of  p  at  low  carrier  densities  are  carried  out 
as  a  function  of  Cd  concentration  (x)  in  Hgi^CdjTe  alloys  and  are  shown  in  Figure  5.  The 
wavelength  dependence  is  similar  in  all  alloys.  However,  as  the  value  of  x  increases,  the 
band  gap  increases  and,  consequently,  p  decrease  with*.  This  is  clearly  seen  in  Figure  5. 

Free  carrier  absorption 

Along  with  TP  A,  further  absorption  light  by  the  photo-generated  carriers  is 
possible.  When  a  hole  is  photo-generated  in  the  heavy-hole  (HH)  VB,  the  electrons  in  the 
light-hole  (LH)  band  can  absorb  one  photon  and  get  excited  to  fill  that  hole.  This  direct 
process  is  normally  known  as  FCA  by  holes.  The  associated  photo-generated  electron  in 
the  CB  also  will  absorb  a  photon.  However  to  satisfy  the  momentum  and  energy 
conservation  condition,  the  photon  absorption  should  accompany  a  phonon  absorption  or 
emission  in  this  indirect  process.  The  free  carrier  absorption  rate  calculated  from  second 
order  perturbation  theory  is  used  to  obtain  FCA  coefficient  ctf  and  FCA  cross  section 
a.  Although  the  basic  formalism  has  been  described  previously,  here  we  provide 
explicit  expressions  for  two  contributions — direct  term  and  indirect  term.  This  helps  us 
separate  the  contributions  from  electrons  and  holes,  for  an  appropriate  use  in  the  rate 
equations.  In  the  direction  term,  electron  from  a  lower  band  absorbs  a  photon  and 
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transfers  to  a  higher  band.  The  initial  and  final  state  wave  vector  remains  unchanged,  but 
the  final  state  energy  is  larger  than  the  initial  state  energy  by  .  We  have, 
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Since  the  energy  separation  between  LH  and  HH  bands  vary  from  0  to  several  electron 
volts,  this  term  is  non  zero  for  hole-initiated  FCA.  However,  the  CBs  are  separated  by 
energies  much  larger  than  the  usual  phonon  energies  of  interest,  and  the  direct  excitation 
of  electrons  from  the  lowest  CB  to  the  next  highest  CB  is  not  possible.  Hence,  in  the 
expression  for  ap,  both  a  and  p  denote  valence  bands  and  the  FD  function  is  evaluated 
with  EfP.  An  additional  mechanism  to  supply  both  momentum  and  energy  is  required  for 
CB  electrons  to  participate  in  the  FCA.  The  most  dominant  contribution  usually  arises 
from  longitudinal-phonon  (LO)  assisted  FCA.  Using  Frohlich  Hamiltonian  to  describe 
electron-phonon  interaction,  the  absorption  coefficient  for  this  indirect  process  can  be 
written  as, 
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where  the  ftcoLO  and  Nlo  are  LO  phonon  energy  and  population,  k0  and  k„  are  zero  and 

infinite  frequency  dielectric  constant,  and  ±  sign  denotes  phonon  absorption  or  emission. 
Other  symbols  have  the  same  meaning  as  those  defined  in  the  previous  section.  For  the 
conduction  electrons,  the  initial,  intermediate  and  final  states  will  all  be  in  the  lowest  CB. 
However  for  the  holes,  the  initial  and  final  states  can  be  in  either  of  LH  or  HH  bands.  It 
should  be  noted  that,  depending  on  the  bands  that  contribute  to  the  absorption,  the  FD 
functions  are  evaluated  using  appropriate  quasi  Fermi  level.  The  total  FCA  coefficient,  af 
is  sum  of  the  two  contributions.  The  FCA  cross  section  for  electrons  (holes),  oe(h)  is  ratio 
of  the  coefficient  to  the  electron  (hole)  density,  Ne(h).  Simply, 
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af  =aD  +alD 


As  expected,  the  calculated  rate  for  FCA  by  electrons,  which  has  only  indirect 
term,  is  much  smaller  than  that  for  the  holes.  In  all  our  calculations,  we  found,  in 
agreement  with  literature,  that  the  FCA  coefficient,  for  both  electrons  and  holes,  is 
linearly  dependent  on  the  carrier  density.  Consequently  ct  is  independent  of  the  carrier 
density.  The  calculated  ct  is  plotted  as  a  function  of  wavelength  for  the  four  materials  in 
Figure  6  for  InAs  and  HgCdTe  alloy.  Note  that  the  electron  contribution  is  much  smaller 
than  the  hole  contribution  as  expected.  The  total  cross  section  is  small  at  shorter 
wavelengths  because  of  the  unavailability  of  holes  at  those  energies  in  the  HH  band. 

Since  the  hole  density  is  more  near  the  zone  center,  the  FCA  cross  section  increases  with 
the  photon  wavelength.  However,  because  of  the  reduced  intermixing  of  symmetries  in 
both  LH  and  HH  band  states  near  the  zone  center,  the  dipole  matrix  element  decreases 
(and  reaches  zero  at  the  center).  Consequently,  the  cross  section  reduces  with  further 
increase  in  wavelength.  This  behavior  is  seen  in  all  four  materials  considered  here,  but 
shown  only  for  InAs  and  HgCdTe  alloy  in  Figure  6. 

When  light  energy  is  absorbed  in  semiconductors,  equal  densities  of  electrons  and 
holes  are  created,  which  in  turn  release  that  energy  to  the  lattice  in  the  form  of  heat.  To 
model  light  propagation  through  semiconductors,  it  is  therefore  necessary  to  know'  the 
FCA  cross  section  as  a  function  of  both  carrier  concentration  and  temperature.  The  ct 
calculated  for  four  materials  as  a  function  of  wavelength  and  temperature  are  shown  in 
Figure  7.  In  the  literature,  ct  is  usually  assumed  to  be  independent  of  temperature  T. 
However,  Figure  7  shows  that  ct  varies  strongly  with  T  —  at  long  wavelength  ct 
decreases  with  T  while  at  short  wavelength  it  increases  with  T.  This  dependence  can  be 
explained  using  the  second-order  perturbation  theory.  The  absorption  coefficient  is 
approximately  the  integrated  product  of  the  hole  density  and  optical  matrix  element. 
Although  this  coefficient  is  divided  by  the  total  hole  density  to  get  ct,  only  a  fraction  of 
the  holes  participate  in  the  absorption  process  because  of  energy  and  momentum 
conservation  conditions.  Even  when  the  matrix  elements  do  not  change  with  T,  the  ratio 
of  participating  holes  to  the  total  number  of  holes  has  strong  T  dependence.  For  larger 


10 


photon  energies  (or  shorter  wavelengths),  the  VB  states  away  from  the  zone  center 
participate  in  absorption  and  the  hole  density  in  those  states  increase  with  T.  For  longer 
wavelength  photon  absorption,  the  holes  near  the  zone  center  participate  and  the  hole 
density  decreases  with  T  as  the  holes  are  thermally  excited  to  higher  energies. 
Consequently,  the  a  decreases  with  T. 

We  have  calculated  the  T-dependence  of  a  in  InAs  and  Hgo.66Cdo.34Te  at  a 
wavelength  of  at  4.8  pm  and  in  InSb  and  Hgo.76Cdo.24Te  at  a  wavelength  of  at  9.6  pm. 
Figure  8  shows  that  the  a  at  9.6  pm  decreases  by  a  factor  of  14  to  16  as  the  T  is  increased 
from  77K.  to  500K.  cr  at  4,8  pm  in  Hgo.66Cdo.34Te  is  nearly  a  constant,  but  increases 
rapidly  with  T  in  InAs.  It  is  important  to  note  that  T  dependence  of  a  is  in  general  very 
strong  except  at  one  or  two  wavelengths  where  the  hole  density  does  not  change  with  T. 

The  wavelength  dependent  a  is  calculated  as  a  function  of  x  in  Hgi^Cd^Te  alloys  at 
T=300K  and  is  shown  in  Figure  9.  At  smaller  wavelengths,  cr  is  independent  of  the  x,  as 
the  participating  hole  population  does  not  vary  much  with  x.  As  x  is  increased,  the  spin- 
orbit  coupling  increases  and  the  LH  band  bends  more  away  from  HH  band."0 
Consequently,  the  energy  conservation  condition  is  satisfied  in  the  BZ  closer  to  zone 
center.  The  resulting  increase  in  participating  hole  density  causes  a  to  increase  with  x  as 
seen  in  Figure  9. 

Lifetimes 

The  carriers  excited  by  the  light  absorption  can  participate  in  further  absorption  of 
light  before  they  recombine.  Hence  the  lifetimes  t  of  the  excited  electrons  play  a  vital 
role  in  light  propagation  through  the  material.  The  full  band  structures  and  accurate 
scattering  matrix  elements  are  used  in  the  lifetime  calculations  as  described 
previously.25,36’37  However,  those  calculations  are  now  generalized  to  include  the  quasi 
Fermi  levels  for  electrons  and  holes.  The  lifetimes  limited  by  the  Auger  and  radiative 
recombination  are  obtained  as  functions  of  temperature  and  carrier  densities. 

Parameter  fit 
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3  6  22 

Although  the  light  propagation  modeling  studies,  '  ’  including  solutions  to  the 
rate  equations,22,24  use  only  a  few  parameters—  p,  a,  x,  N  (-no+An),  T,  and  thickness  d, 
the  complication  arises  because  of  their  interdependence.  For  example,  p  is  a  function  of 
both  N  and  T,  but  N  depends  on  p  and  T.  Similarly  t  depends  on  T  and  N,  which  in  turn 
depends  on  T.  Eventually,  p,  a  and  x  will  determine  the  increase  in  T,  which  will  in  turn 
decide  the  values  of  those  input  parameters.  Hence,  for  a  self-consistent  solution  to  the 
rate  equations,  it  is  preferable  to  have  a  closed  functional  form  for  the  parameters  p,  cr 
and  x  which  is  valid  in  the  temperature  and  carrier  density  ranges  of  interest.  The  values 
of  p  and  a  have  been  calculated  at  a  wavelength  of  4.8  pm  in  InAs  and  Hgo.66Cdo  34Te 
and  at  9.6  pm  in  InSb  and  Hgo.76Cdo.24Te,  for  77K  <  T  <  500K  and  1014  cm°  <  An  <  101S 
cm'3.  We  have  them  fitted  the  values  to  the  following  expressions. 


The  fitted  parameters  are  given  in  Table  1  for  four  materials  considered  here.  The  table 
lists  parameters  only  for  AR  lifetimes.  We  had  previously  calculated  radiative 
recombination  lifetimes  and  found  them  in  general  to  be  very  long  for  these  materials.^ 
In  addition,  because  of  photon  recycling  present  in  thick  samples  studied  here,  the  RR 
lifetimes  are  extremely  long  to  have  negligible  effect  on  light  propagation  and  hence  the 
parameters  are  not  listed  here.  The  AR  lifetimes  in  Hgo.76Cdo.24Te  are  not  calculated  and 
not  listed,  as  they  are  not  expect  to  be  substantially  different  from  that  in  Hgo.66Cdo.34Te. 


EXPERIMENTS  AND  RESULTS 

In  the  approximation  that  the  change  in  refractive  index  caused  by  photoexcited 
carriers  in  the  material  is  negligible,  the  Helmholtz  equation  can  be  reduced  to  a  simple 
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rate  equation.  A  complete  and  self-consistent  solution  to  these  rate  equations  is  needed 
to  develop  a  detailed  understanding  of  the  mechanisms  affecting  the  light  propagation  in 
the  material.  Recently  we  used  above  calculated  parameters  for  InAs  and  solved  the  rate 
equations  to  explain  the  results  obtained  in  our  pump-probe  experiment.38  Since  the 
details  are  published  else  where/ 8  only  the  results  are  discussed  briefly  here  to  provide 
continuity. 

in  the  pump-probe  experiment  the  charge  carriers  were  generated  by  the  pump 
using  TPA.  These  carriers  absorb  the  probe  photos,  resulting  in  very  small  probe 
transmission.  As  the  carriers  recombine,  the  probe  transmission  increases  and  reaches  its 
incident  value.  The  TPA  used  here  ensures  that  the  entire  thickness  of  the  sample  is 
utilized  for  carrier  generation  at  low  incident  irradiance.  The  4.8  pm  pump  source  was 
obtained  by  frequency  doubling  a  TEA  CO2  laser  to  provide  128  ns  (full-width  at  e  'of 
the  maximum)  duration  pulses  with  energies  up  to  10  mJ.  The  measurements  are  carried 
out  at  room  temperature.  First,  we  carried  out  a  set  of  nonlinear  absorption  experiments 
with  only  the  pump  present.  The  measured  transmitted  intensity  as  a  function  of  incident 
intensity  and  sample  thickness  is  shown  in  Figure  10.  As  the  intensity  is  increased,  more 
carriers  are  generated  and  the  absorption  increases,  resulting  in  a  decreased  transmission. 
Similarly,  as  the  thickness  is  increased,  light  is  absorbed  more.  These  trends  are  clearly 
observed.  Also  in  Figure  10,  we  have  plotted  the  calculated  transmitted  intensity  obtained 
from  solving  the  rate  equations  with  the  values  given  in  Fable  1 .  In  addition  to  the  Auger 
recombination,  we  have  included  lifetimes  limited  by  Shockley-Read-Hall  (SRH) 
mechanism.  The  SRH  scattering  is  not  intrinsic  to  the  material  and  depends  on  external 
variables  such  as  the  growth  temperature  and  pressure.  Our  detailed  calculations  carried 
out  previously39  indicate  a  nearly  T-independent  SRH  liletime  of  200  ns  in  a  2  x1016  cm 
n-doped  InAs.  Although  the  InAs  sample  studied  here  is  only  1015  cm'3  n-doped,  the 
excess  carriers  generated  by  the  pump  is  of  the  order  ot  1-2  x  1016  cm  3,  and  justifies  the 
use  of  the  above  calculated  SRH  lifetimes.  We  further  assumed  that  SRH  lifetime 
decreases  linearly  with  excess  carrier  density.  We  see  that  calculated  transmission 
intensity  agrees  very  well  with  measured  values  for  all  thicknesses  studied.  The 
calculations  are  then  extended  to  understand  the  pump- probe  results  shown  in  Figure  1 1 . 
With  only  the  SRH  values  used  in  an  approximation,  we  see  that  the  time-dependent 
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transmitted  probe  intensity  calculated  for  the  three  input  energies  explain  the  observed 
variation  very  well. 

CONCLUSIONS 

We  have  used  a  full  bandstructure,  obtained  from  a  hybrid  pseudopotential  tight- 
binding  Hamiltonian,  to  calculate  the  temperature  and  carrier  density  dependent  two 
photon  absorption  coefficients,  free  carrier  absorption  cross  section,  and  Auger  lifetimes 
in  InAs,  InSb  and  HgCd'l'e  alloys.  Accuracy  of  the  calculations  are  enhanced  with  use  of 
the  quasi  Fermi  levels,  intrinsic  carrier  densities,  detailed  optical  and  scattering  matrix 
elements,  and  complete  integration  over  the  entire  BZ.  The  calculated  parameters  are 
fitted  to  closed-form  functions  and  used  in  the  rate  equations  to  obtain  the  transmitted 
intensity  as  a  function  of  time,  incident  intensity,  and  sample  thickness.  The  transmission 
spectra  calculated,  with  only  SRH  lifetimes  as  a  parameter,  agree  very  well  the  measured 
values.  We  showed  that  with  accurate  evaluation  of  all  parameters,  the  modeling  can 
provide  a  detailed  guidance  in  understanding  the  underlying  mechanisms  and  in  choosing 
materials  for  light  propagation  applications.  The  closed-form  expressions  for  the 
fundamental  parameters  given  here  can  be  used  in  other  accurate  light  propagation 
models. 
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TABLE  1 


Parameters 

InAs 

Hgo.6eCdo.34Te 

InSb 

Hgo,76Cdo,44Te 

ao 

0.20516 

2.1620 

2.5913 

48.074 

ai 

0.0015299 

-0.00394730 

0.0039691 

-0.45516 

a2 

-1.859x1 0"6 

1.3372xl0‘5 

0.00012589 

0.0029993 

a3 

1.8466xl0‘9 

-3.1219xl0’8 

-4.6698X10'7 

-8.3954xl0'6 

a4 

3.4373x1  O'12 

3. 3459x1  O'" 

6.6038x1  O*!0 

7.71 8x1  O'9 

bo 

-0.034375 

0.11742 

-0.0019728 

0.3932 

b] 

0.0020645 

0.00041377 

0.00031268 

-0.0033771 

b2 

-2.1 108x1 0*6 

5.3520xl0'6 

5.2291  xlO'7 

2.0397x1 0"5 

b3 

-6.5561xlO'10 

-6.65 17x1  O’9 

7.3333X10'10 

-5.051x10-* 

b4 

5.6352xl0-'2 

5.4334xl0‘12 

3.0668x1  O'12 

4.563x10*" 

co 

0.45859 

1.2756 

-0.04522 

0.21089 

Cl 

-0.0013534 

-0.00095861 

0.0 

0.0032988 

C2 

1.0666xl0’s 

1.0237xl0"5 

0.0 

-3.8836X10'5 

C3 

-1.1 906x10'* 

-1,8417x10-* 

0.0 

1 . 1 845x  1 0"7 

C4 

4.3216xl0‘12 

-7.1237X10'12 

0.0 

-1.1285xlO*10 

CF0 

-2.8164 

-0.95501 

46.1 

65.296 

Oi 

0.044098 

0.16677 

-0.2929 

-0.57162 

ct2 

a3 

o4 

5.9962e-06 

-2.2463e-07 

-0.00081382 

1.4824X10’6 

0.00082667 
- 1 . 1092x  1  O'6 

0.0021109 

-3.4952xlO’6 

2.3863e-10 

-9.422X10'10 

5.6047xl0'10 

2.1284xl0'9 

go 

40.467 

43.429 

47.74 

43.429 

gl 

-0.02737 

-0.06712 

-0.15205 

-0.06712 

g2 

-8.597x1  O'6 

2.480x1  O’4 

0,0007528 

2.480x1 0"4 

g3 

2.2386x10‘7 

-3.7801  xlO-7 

-1.6296x1 0‘6 

-3.7801x1  O'7 

g4 

-3.122X10'10 

1.8071X10'10 

1 .241  xlO"9 

1.8071x10"'° 

ho 

14.336 

15.1230 

14.26 

15.1230 

hi 

0.037678 

0.034190 

0.038101 

0.034190 

h2 

-0.00011995 

-1.1134xl0'4 

-0.00010939 

-1.1134xl0"4 

h3 

1.8793e-07 

1.8605X10’7 

1.4327xl0"7 

1.8605xl0'7 

h4 

-1.152e-10 

-1.2588x10-'° 

-6.6935x10-“ 

-1.2588x10-'° 

Po 

9.2435 

11.534 

8.5016 

8.7934 

Pi 

-2.7345 

-3.9967 

-1.7296 

-2.1725 

P: 

0.0 

0.0 

0.0 

0.19217 

P3 

0.0 

0.0 

0.0 

-0.013391 

P4 

0.0 

0.0 

0.0 

0.00036021 
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Figure  Captions: 

Figure  1.  Variation  of  energy  gap  (dashed  line)  and  intrinsic  Fermi  energy  (dotted  line) 
with  temperature.  Also  shown  are  electron  quasi  Fermi  energy  (thin  solid  line),  hole  quasi 
Fermi  energy  (thick  solid  line)  for  various  photoexcited  carrier  densities  (An=Ap). 

Figure  2.  Variation  of  intrinsic  carrier  density  with  temperature 

Figure  3,  Two  photon  absorption  coefficient  as  a  function  of  wavelength  and 
temperature. 

Figure  4,  Two  photon  absorption  coefficient  as  a  function  of  excited  carrier  density  at  (a) 
at  300K  and  (b)  77K. 

Figure  5.  Two  photon  absorption  coefficient  as  a  function  of  wavelength  in  HgCdTe 
alloys  at  300K. 

Figure  6,  Wavelength-dependent  free  carrier  absorption  cross  section  at  300K 
Figure  7.  Wavelength-dependent  free  carrier  absorption  cross  section  at  77K  and  300K. 

Figure  8.  Temperature-dependent  free  carrier  absorption  at  4.8  pm  in  InAs  and 
Hgo  66Cdo.34Te  and  at  9.6  pm  in  InSb  and  tlgo.76Cdo.24Te  at  300K, 

Figure  9.  Wavelength-dependent  free  carrier  absorption  cross  section  as  a  function  of  Cd 
concentration  in  HgCdTe  alloys  at  300K.. 

Figure  10.  Comparison  of  normalized  pump  transmission  coefficient  measured  (data)  at 
300K.  in  1015  n-doped  InAs  with  that  calculated  (solid  lines)  for  four  thicknesses. 

Figure  1 1.  Comparison  of  normalized  probe  transmission  coefficient  measured  (data)  at 
300K  in  1015  n-doped  InAs  with  that  calculated  (solid  lines)  for  three  beam  energies. 
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Figure  1  Krishna  mu  rthy  et  al. 
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Figure  2  Krishna  mu  rthy  et  al 
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Figure  3  Krishnamurthy  et  al. 
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Figure  4  Krishnamurthy  et  at 
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Figure  5  Krishnamurthy  et  al. 
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Figure  6  Krishnamurthy  et  al. 
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Figure  7  Krishnamurthy  et  al 
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Figure  8  Krishna  mu  rthy  et  al. 
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Figure  9  Krishnamurthy  et  aL 
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Figure  10  Krishnamurthy  et  al. 
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Figure  11  Krishna  mu  rthy  et  al. 
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